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Abstract. We report on the application of a cascade -I- viscous hydro -|- cascade model 
for heavy ion collisions in the RHIC Beam Energy Scan range, ^snn = 6.3... 200 GeV. By 
constraining model parameters to reproduce the data we hnd that the effective(average) value 
of the shear viscosity over entropy density ratio rj/s decreases from 0.2 to 0.08 when collision 
energy grows from ^snn « 7 to 39 GeV. 


1. Introduction 

The principal goal of the Beam Energy Scan program at Relativistic Heavy Ion Collider (RHIC) 
is to study the transition from baryon-free to baryon-rich regime of strongly interacting QCD 
matter by decreasing the collision energy. One of questions is what the transport properties of 
the medium with high baryon density, created at lower beam energies, are. 

To address this question we use a state-of-the-art event-by-event viscous hydrodynamic -|- 
cascade (or hybrid) model to describe the evolution of the collision system. 

2. The model, its parameter space and results 

The hybrid model consists of the pre-equilibrium stage, the fluid dynamical stage and the 
hnal hadronic cascade. The pre-thermal stage of evolution is described by UrQMD cascade 
mm- The interactions are calculated until a hypersurface of a constant Bjorken proper time 
T = Vf^ - z‘^ = To, which is a parameter of the model. At r = tq the energy, momentum and 
charges of particles are smoothly distributed to fluid cells according to a Gaussian prohle. The 
contribution of each particle to a fluid cell {i,j, k} is given by: 

= {P^V»} -C-exp ■ (1) 

where Ax*, Ayj, Arjk are the differences between particle’s position and the coordinates of the 
hydrodynamic cell {i,j, k}. Due to the lack of knowledge about thermalization process, R± and 
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Figure 1. Parameter dependence of the total yield at midrapidity (left) and pT integrated 
elliptic flow (right) in 20-30% central Au-Au collisions at y^SNN = 19-6 GeV. Figures are from 

Ref. [6]. 


Rri are another free parameters of the model. These distributions are used as an initial state 
for 3-dimensional hydrodynamic expansion, described by relativistic viscous hydrodynamic code 
vHLLE [3j. We use the equation of state based on Chiral model |lj in the fluid stage. We set 
bulk viscosity and baryon charge diffusion to zero in the fluid phase, and study the effects of 
shear viscosity only. When the energy density reaches a certain value e = Cgw) the hydrodynamic 
medium is converted to particles (particlized) according to Cooper-Frye prescription [5] . The last 
stage of evolution - hadron rescatterings and decays - is described again with UrQMD cascade. 
For more details see Ref. [6j. 

As a result, the model has three parameters related to the initial state (starting time tq and 
Gaussian smearing radii i?_L, i?r?), one parameter of the hydrodynamic medium (shear viscosity 
over entropy ratio r]/s) and one parameter related to switching from hydrodynamics to cascade 
(switching energy density Cgw)- The main goal is to extract the value of the physical parameter 
r]/s from comparison with available experimental data. This goal is, however, complicated by 
the fact that the observables in the model depend significantly on all the mentioned parameters, 
not only rj/s. A somewhat stronger sensitivity to the initial state parameters compared to 
hydrodynamic studies at full RHIG or LHG energies [someone] is explained by shorter duration 
of hydrodynamic phase at lower beam energies. 

To learn how exactly final observables depend on the model parameters, we consider the setup 
with R± = Rrj = 1.0 fm, tq = 2R/{'yVz), p/s = 0 and Cs* = 0.5 GeV/fm^ as a “default” scenario. 
On Fig. we show how deviations of parameter values from the “default” scenario affect the 
multiplicity of all hadrons at midrapidity and pT integrated elliptic flow of all charged hadrons 
at one particular collision energy .^snn = 39 GeV. One can summarize the dependencies seen 
in Fig. [^in a following way: 

• Increased R± smoothens the initial energy density profile in the transverse plane, which 
leads to smaller gradients and less explosive transverse expansion. Increased R± also 
results in decreased ellipticity an initial energy density profile, which is hydrodynamically 
translated into smaller final elliptic flow. 

• In a similar manner, the increase of Rrj leads to shallower longitudinal gradients and 
weaker longitudinal expansion. Thus more energy remains at midrapidity to form stronger 
transverse expansion, which increases V 2 - On the other hand, larger also results in larger 
initial entropy of the fluid, which considerably increases the final particle multiplicity. 

• Increased tq leads to a shorter lifetime of the hydrodynamic phase, as a resnlt of longer 
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Table 1. Collision energy 

dependence of the model parameters, 
which allows the model to reproduce 
the experimental data in BES region 
and higher RHIC energies. 


pre-thermal phase. At the same time tq enters the Gaussian energy/momentum smearing 
profile. Thus its increase acts opposite to the increase of The results show that the 
second effect is stronger. 

• Increased Cgw shortens the effective lifetime of the hydrodynamic phase. The shorter time to 
develop radial and elliptic flows is not fully compensated by the longer cascade phase, which 
results in smaller final V 2 - Since the total entropy is conserved in the ideal hydrodynamic 
expansion, but increases in the cascade stage, the final particle multiplicity increases with 
the increase of Cgw 

The found dependencies allow us to adjust the values of parameters individually for every 
collision energy to approximately reproduce the available experimental data points in the BES 
region. To simplify the analysis, we keep the same value of switching energy density Csw at all 
collision energies. The adjusted values of the parameters are shown in Table Some of the 
results are shown in Eig. 

The shear viscosity over entropy density ratio 77 /s in fluid phase is taken to be independent 
of temperature and chemical potentials, and changes only with collision energy. Therefore the 
value of ry/s can be thought of as effective one for given computational setup. The parameter 
adjustment procedure yields 77/5 values which decrease with increasing collision energy. At lower 
collision energies the fluid phase probes thermodynamic region of lower temperature, limited 
by the particlization temperature from below, and higher chemical potential. One expects 
that physical 77/5 depends on temperature and chemical potential, and not on collision energy. 
Therefore we conclude that the .^snn dependence of the effective rj/s we found hints that in more 
realistic simulations {rj/s){T, hb) should increase with increasing baryon chemical potential ^b- 

3. Summary 

To summarize, the 3 dimensional cascade + viscous hydro + cascade model is applied for the 
description of heavy ion collisions in the collision energy range of RHIC Beam Energy Scan 
program, .^/snn = 6.3 ... 200 GeV. The parameter space of the model is explored. The values of 
model parameters are adjusted for different collision energies in order to approach the existing 
experimental data for rapidity and px distributions, and elliptic flow coefficients of hadrons. 
As a result it is found that the value of effective shear viscosity over entropy density ratio r]/s 
decreases from 0.2 to 0.08 when the center of mass collision energy increases from 7.7 to 39 GeV, 
and stays at 0.08 for .^snn > 39 GeV. 
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Figure 2. Rapidity (top left) and transverse momentum (bottom left) distributions of pions 
and kaons at Fiab = 40 A GeV Pb-Pb collisions, pseudorapidity distributions of all charged 
hadrons at ^/s = 19.6,62.4 and 200 GeV Au-Au collisions (top right), and px integrated flow 
coefficients for different collision energies [9] (bottom right), all calculated in the model and 
compared to the experimental data from NA49 [7] and PHOBOS [8]. Figures are from Ref. [ 6 ]. 
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